Anisotropic stress and stability in modified gravity models 
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The existence of anisotropic stress of a purely geometrical origin seems to be a characteristic of 
higher order gravity models, and has been suggested as a probe to test these models observationally, 
for example in weak lensing experiments. In this paper, we seek to find a class of higher order gravity 
models of f{R, G) type that would give us a zero anisotropic stress and study the consequences for 
the viability of the actual model. For the special case of a de Sitter background, we identify a 
subclass of models with the desired property. We also find a direct link between anisotropic stress 
and the stability of the model as well as the presence of extra degrees of freedom, which seems to be 
a general feature of higher order gravity models. Particularly, setting the anisotropic stress equal 
to zero for a de Sitter background leads to a singularity that makes it impossible to reach the de 
Sitter evolution. 

PACS numbers: 04.50.Kd, 98.80.-k, 95.36.+X 



I. INTRODUCTION 



Even though it is now well over ten years since super- 
nova data [ir 01 led to the general acceptance amongst 
cosmologists that the expansion of the Universe is accel- 
erating, there are no natural models that could explain 
this phenomenon. The nature of the postulated dark en- 
ergy which is responsible for the accelerated expansion 
is still a mystery ^-i9j- This fact, as well as expecta- 
tions from candidates for the "theory of everything" , like 
string theory, have led researchers to investigate mod- 
els in which General Relativity (GR) itself is modified 
One crucial question therefore is whether it is 
possible to distinguish between a dark energy "field" or 
fluid on the one hand, and an accelerated expansion due 
to a modification of the theory of gravity on the other 
hand. 

Although strictly spea king cosmological probes cannot 
provide conclusive proof [19l - l2l| . there are certain signa- 
tures that can point the way. One of them is the presence 
of a significant anisotropic stress: canonical scalar fields 
do not create additional anisotropic stress, while the 
modified-gravity (MG) models like scalar-tensor theories, 
brane-world models like the Dvali-Gabadadze-Porrati 
(DGP) model [22] and f{R, G) type theories generically 
induce a large effective anisotropic stress. 

In this paper we investigate one specific class of models, 
f{R^ G) type modifications of GR, and ask the question 
whether it is possible to construct viable models with a 
vanishing, or arbitrarily small effective anisotropic stress. 
Or in other words, is it possible to mimic "GR" with these 
models, at least up to first order in perturbation theory 
and in the sense that the extra anisotropic stress is small 
enough? Since f{R) models have many things in common 



with scalar-tensor theories, we expect that our discussion 
is also relevant for those models, and as we discuss later, 
also for DGP and other braneworld models. 

We structure the paper as follows: The following sec- 
tion serves to define our notation and introduces the 
main equations. In Section Hill we discuss the notion of 
anisotropic stress in general, and how this plays an im- 
portant role in modified gravity models and then we in- 
vestigate the possibility of a vanishing anisotropic stress 
in the particular cases of f{R) and /(G) models, before 
we look at the more general f{R,G) case. In Section 
IIVI we identify and discuss the link between anisotropic 
stress and stability in modified gravity models in the con- 
text of both homogeneous and inhomogeneous perturba- 
tions around de Sitter space. We further derive the rel- 
evant stability conditions. We generalize the discussion 
to arbitrary backgrounds in Section |V] and give some re- 
sults for a matter dominated evolution. In Section IVll we 
apply the above to characteristic toy models, and then 
discuss our conclusions. Some explicit intermediate cal- 
culations and formulas can be found in the Appendices. 



II. ACTION AND EQUATIONS OF MOTION 

We are interested in the class of models described by 
the action 



S = 



d^xy/^[f{R,G) +L 

matter] ; 



(1) 



M 
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where R and G are the Ricci and Gauss-Bonnet scalar 
respectively, M denotes the four dimensional spacetime, 
and Lmattcr is the Lagrangian for any matter fields or 
fiuids present. "/" is an analytic function and its form 
is constrained by both classical and quantum stability 
requirements as well as agreement with both large scale 
and solar system data. (We will revisit this point in Sec- 
tion |lVl) In the following, we will work in natural units 
where SttG = c = 1, unless otherwise stated. 
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The Gauss-Bonnet term, which is a topological invari- 
ant in four dimensions, is defined as 



(2) 



and belongs to the so-called Lovelock scalars [23], which 
is the first term in a series of non-linear corrections to 
the Einstein-Hilbert term. 

The class of models described by action ([T|) are a nat- 
ural generalisation of the well-studied J{R) models [lol - 
ITsj . The introduction of curvature scalars other than the 
Ricci one has been physically motivated by the lower en- 
ergy limit of string theory 24], although care must be 
taken to avoid generic instabilities (e.g. [IH). It is well 
known that actions like ([l} can be expressed in an equiv- 
alent way as a scalar-tensor action in the so-called Jordan 
frame, however, in this paper we will work in the original 
representation. 

Varying action ([Ij with respect to the metric g^'', and 
usin g th e Bianchi identities, we get the equations of mo- 
tion [H, [131 



Ffn —In Rl — rTi(mattor) , rp{cS) 



(3) 



where TjfJ^"^ is an effective energy-momentum tensor of 
geometrical origin (in contrast to the ordinary matter 
one) and is defined as 

Tjf^ = - (V^V.F - g^.DF + 2i?V^V.e 



where we used the additional definitions 

_ df{R,G) 
^-f^= OR ' 



df{R,G) 
dG '■ 



(4) 



(5) 



(6) 



V{R, G)=RF + iG- f{R, G). 



Taking the limits <^ -> and F ^ 1 in ([3]), we recover 
the /(i?) and R + f{G) equations of motion respectively. 

We can formally recover a GR-like equation of mo- 
tion by moving additionally {F — l){R^i, — ^g^^R) to the 

right hand side of Eq. (jS]) and adding this to tI^^J^^ to 
form a total effective "dark energy" energy-momentum 
tensor, j^^J^^^"^^^) _ ca,n then compute effective fluid 
quantities, for example the equation of state, pressure 



perturbation and anisotropic stress, that a dark energy 
needs to have to lead to the same cosmological observa- 
tions as the original /(i?, G) model. But strictly speak- 
ing, although covariantly conserved, is not an 
energy-momentum tensor in the usual sense, since it is 
a function of the spacetime geometry and its first and 
second derivatives. 

In the following, we will be interested in homogeneous, 
isotropic and flat cosmologies, described by the fiat, four 
dimensional FRW metric 



(8) 



with a{t) the scale factor. In this background, the two 
key quantities, R and G, can be expressed purely as a 
function of the Hubble parameter H = H{t) and its time 
derivative, 



6 2H^ + H 



G{t) = 2AH^ [H 



H 



(9) 
(10) 



The t-t component of the f{R,G) equations of motion 
([3]) gives a modified version of the usual Friedman equa- 
tion which reads as 

with dots denoting differentiation with respect to cosmic 
time t. Notice that the above equation is of fourth order 
with respect to the scale factor, in contrast to the usual 
Friedman equation. 



III. THE EFFECTIVE ANISOTROPIC STRESS 
IN HIGHER ORDER GRAVITY 

Let us here introduce the notion of anisotropic stress 
in gravity. As a starting point, we consider scalar pertur- 
bations around a flat FRW background in the conformal 
Newtonian gauge, where the metric is of the forni^ 



ds' 



-(1 + 2il;)dt^ + a{tf (1 - 2(1)) dx^ 



(12) 



and the gravitational potentials = V'(x, i) and (j) = 
(7) (/)(x, t) are closely related to observations: light deflec- 
tion is sourced by the lensing potential (f) + ip and non- 
relativistic particle motion by alone. 

The scalar anisotropic stress H is then defined as the 
difference in the potentials 



V' = n(x,t), 



(13) 



or the difference of the relevant potentials in some other 
gauge. Equation (jl3p is called the anisotropy equation, 



^ Here we follow the notation of [27|| . 



^ The general form of the perturbed line element is given in Ap- 
pendix |B] 
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and can be found by calculating the ij {i ^ j) component 
of the perturbed equations of motion around the FRW 
metric, 

SG\ - IgiSG''^, = J'7^(cff,total)» . _ 1 »^j,(off,total)K^ 



jj(cff)i 



(14) 



from which one then extracts the scalar part as usual to 
get 



(15) 



We emphasize that this is the anisotropic stress one 
would infer by assuming GR to hold, not only the 
anisotropic stress from the matter fields. Indeed, here 
we are precisely interested in the contribution to (j) — 
due to a modification of gravity. While relativistic par- 
ticles do induce an anisotropic stress, it is small at late 
times and we will neglect the contribution of 
equation ([3]) to — -0. Notice that because of the na- 
ture of in modified gravity theories, the r.h.s of 
above equation will in principle have a spacetime depen- 
dence, i.e it will be a function of 0, ip as well as their first 
and second derivatives with respect to time (in Fourier 
space), in contrast to GR, where the r.h.s is just a func- 
tion of the matter content. The usefulness of (|15l) is that 
it has a GR-like l.h.s., allowing to compute predictions 
for cosmological observations as usual, while all the extra 
contributions are moved to the r.h.s. and interpreted as 
a "modified gravity energy-momentum tensor" . 

In particular, for GR (and neglecting any relativistic 
species) we have n^''^' = and therefore (j) = ip at all 
times. Therefore, the inequality of the Newtonian po- 
tentials is a "signature" of departures from GR on large 
scales ^9j. The ratio <f>/ip, or variables derived from it, 
like ri(t,k) = ^ — 1, can be extracted observationally 
by combining weak lensing experiments with e.g. galaxy 
surveys or redshift space distortions, making cosmologi- 
cal observations a powerful test of GR ^] . Current lim- 
its are rather weak, with deviations of order unity from 
77 = still allowed, but future probes will measure the 
ratio ■0/0 with an accuracy of a few percent (e.g. [29 33] ) . 

In this paper, we raise and investigate the following 
question: Can we construct a viable modified gravity 
model with 0/0 — 1, or in other words, is ^ an 
unavoidable consequence of modifying gravity to explain 
the dark energy? We will try to answer this question 
step by step, by investigating the anisotropy equations 
of f{R), R + f{G) as well as of the more general f{R, G) 
gravity models. 

The equations for general spaces tend to be complex 
and in general do not admit simple solutions. For this 
reason in this paper we will first focus on the case of a 
de Sitter background. On the one hand, solutions that 
explain the observed accelerated expansion usually tend 
towards a de Sitter fixed point, and also the observed 
background expansion requires no deviation from p = —p 
for the inferred dark energy component. On the other 



hand, the equations simplify significantly in this limit, 
which allows us to give explicit solutions that we can then 
discuss in detail. We comment on the behavior for other 
backgrounds in section |Vl but leave a fully general study 
for future work. We nonetheless expect our conclusions 
to be quite generic for models that try to explain the 
dark energy. 



A. The anisotropic stress in f{R) models 

Let us begin with the special case of f(R) gravity, de- 
scribed by the action 



S = 



(16) 



which corresponds to the limit of ^ — > of the general 
f{R, G) models. It is well known that these models are 
characterized by an extra, dynamical scalar degree of 
freedom F, which is proportional to the first derivative 
of /(i?), F = fii{R) = f'{R). This can be exphcitly 
seen by taking the trace of the corresponding equations 
of motion and arriving at a Klein-Gordon type equation 
for F. However, unless the theory is written in the so- 
called Jordan frame, the latter degree of freedom is still 
of geometrical origin. 

The expression for the anisotropic stress for the general 
case of f{R,G) gravity is given in equation For 
f{R) gravity, in the Newtonian gauge, it takes the form 



0- 0> 



SF 



n 



(cff) 



(17) 



which holds for any spacetime, not just de Sitter. 

Since SF — fnii{R)6R, the stress contribution is pro- 
portional to the derivative of the extra scalar degree of 
freedom with respect to R, that is, it depends on the 
evolution of the scalar F = f'{R). Seeking a form for 
the function f{R) that would make II'^'^'^^ vanish at all 
times corresponds to solving the equation fnpi = with 
a general solution f{R) = i? + A, i.e. of all f{R) models 
it is precisely GR that satisfies this equation. In other 
words, the requirement of zero anisotropic stress in f{R) 
theories is equivalent to suppressing the extra degree of 
freedom of the theory, leading to the GR limit. (In 
the Parametrized Post-Friedmannian (PPF) framework 
of Ref. [21], /flfl — >• corresponds to i? — >■ 0, B being a 
parameter introduced to quantify the modification from 
GR). 

Although it is not possible to make n*^"*^) exactly zero 
at all times without reverting back to GR, one can try 
to make it sufficiently small for a given cosmological 
period, by an appropriate choice of the model param- 
eters. This corresponds to setting f^R sufficiently smaU 
for some particular initial conditions and ensuring that 
it stays small, by an approriate choice of model. This 
has been done for example in Ref. 3J|. The price one 



pays is a rapid oscillatory behavior for both the gravita- 
tional potentials and the curvature perturbation. What 
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is more, the amplitude of the latter can grow arbitrarily 
as fnn — >■ 0. We will come back to this this later, when 
we will study the relevant stability conditions and will see 
that this is a general feature of f{R, G) and other mod- 
ified gravity models: the existence of anisotropic stress 
is related to the extra scalar degree of freedom of these 
models, and an attempt to suppress it causes unstable be- 
havior. In the f{R) case, suppression of the extra scalar 
corresponds to fun — > 0. 

As can be seen from equation (|17p , another way to force 

would be to impose the condition 5R = 0. The 



(cff) _ 
R 



n 

crucial difference between 5R = and Jrb. = 0, is that 
the latter is a background requirement, i.e a requirement 
on the particular form of the f{R) action. On the other 
hand, the condition 5R — imposes a dynamical con- 
dition on the potentials 0, if) and their first and second 
time derivatives. If we also take into account that in that 
case the l.h.s implies </> = "0, we find the equation 



IP 
H 



+ 2U = 0, (18) 



which not only in general is unstable, but also fixes the 
perturbation evolution needed to keep SR = 0, which is 
in general incompatible with the desired evolution of the 
universe, e.g. structure formation. In other words, the 
requirement SR — imposes an evolution that in gen- 
eral we do not want. For this reason, what we seek in 
this paper is a condition of the first kind, i.e. a condi- 
tion on model space rather than on the evolution of the 
perturbations. 



B. The anisotropic stress in /(G) models 

Since f{R) models do not allow for a vanishing 
anisotropic stress, we will instead look at the other lim- 
iting case of /(i?, G) models, namely those described by 
the action 



S = 



/(G)] 



(19) 



These models posses an instability in the presence of 
a matter fluid, irrespective of the form of the function 
/(G) [sBj, which rules them out as realistic scenarios, 
but here we just want to see whether it is possible to 
construct f{G) models that contribute no additional ef- 
fective anisotropic stress. 

The first term in the action does not contribute any 
extra anisotropic stress. In an FRW background, these 
models posses an extra scalar degree of freedom, propor- 
tional to ^ = /g(G) = /'(G). The anisotropy equation 
in a general spacetime in this case reads as 



(off) 



-I- ij) 



S^, (20) 



with 5^ = fcG^G. For a de Sitter background, the equa- 
tion simplifies to cj) — ip — AHq/qg^G. One possibility 



to have no anisotropic stress is to set fee = at all 
times, leading to the model /(G) = G + A. In four di- 
mensions G is a topological invariant [l^l, i.e. it is a total 
derivative and so it has no contribution to the equations 
of motion, and we are left only with i? -(- A for the rele- 
vant gravitational Lagrangian, which is equivalent to GR. 
Alternatively we require SG — which suffers from the 
same problems as dR = and does not allow in general 
for a sensible evolution of the perturbations. 

For a general background, the similarity to the case of 
f{R) is spoiled by the first two terms in the anisotropy 
equation. In general, the condition on the evolution of (f> 
and 'ip imposed by those terms will again be difficult to 
enforce as a function of time. On the other hand, if the 
background quantities vary only slowly, ^, ^ « 0, then 
the anisotropy equation can be simplified as 



where we used the relation 



H 3 ,^ 

Jp ~ -2(1 +^eff), 



(21) 



(22) 



with WcS = p/p being the effective equation of state pa- 
rameter for the background evolution. Now, the situation 
is again similar to the one encountered for f{R)'. one has 
either to require either Jqq = 0, SG = 0, or WeS = —1/3. 
As discussed above, the first condition leads to GR (in 
which case automatically ^ = ^ = at all times), while 
the second does not allow for an acceptable evolution 
of the perturbations. The third condition, which corre- 
sponds to the evolution of a universe dominated by curva- 
ture, is also not very relevant given current observational 
results in cosmology. 



C. The anisotropic stress in f{R,G) models 

We saw in the previous sections that the vanishing of 
the anisotropic stress in f{R) and /(G) models corre- 
sponds to either trivial or unphysical situations. We now 
turn to study the more general case of /(i?, G) models. 
Here, the function f{R, G) has two contributions, coming 
from the R- and G- part respectively, and the anisotropy 
equation reads us 



^ F 



5F + 4iJ^V - 4^0 -I- 4 + iJ (5^ 



(23) 



Unlike the f{R) case, where we simply had to demand 
that fB.R.{R) = 0, the nature of the anisotropy equation 
here does again not allow us to write down an explicit 



^ We obtain efFectively the same condition on scales that are well 
inside the horizon, k S> a,H, as (5^ is in general boosted by factors 
of {k/ (aH))'^ relative to if> and -0. 



5 



condition for the function f{R, G) that would give a zero 
anisotropic stress contribution in a general spacetime: as 
in /(G) models, we find extra factors of (/), ip and their 
time derivatives. The only case for which we can find a 
simple condition is for the de Sitter spacetime, and there- 
fore we shall restrict ourselves in this case for the time 
being. Furthermore, for models that try to explain the 
dark energy, it is at late times that we expect modifica- 
tions of gravity to become important, and that deviations 
from GR should appear in observations. For such a late- 
time accelerating epoch, a de Sitter spacetime is expected 
to provide a reasonable approximation. 

The anisotropy equation in de Sitter space reads as 

.nr^ + nr-n£f), (24) 

where, as before, we have defined the contribution coming 
from the R- and G- part of the action respectively as 

n^.f and n^f^^mfj. (25) 

Notice that this case is just the sum of the corresponding 
limiting cases of f{R) and R + f{G) gravity respectively, 
although now either term depends on both R and G. 

We now ask the same question as before: Is it possible 
in this case to find a class of f{R, G) models that give 
a zero anisotropic stress nj^f'' = 0, having at the same 
time a sensible evolution of the perturbations? By in- 
spection of (pSj) one can see that in order for the total 
scalar anisotropic stress to be zero, we require that at all 
times 

I[^f^=-I[^f\ (26) 

In other words, we require that the particular anisotropic 
stress contributions have equal magnitude and opposite 
sign at all times, or at least for the cosmological era of 
interest. 

We can rewrite condition using the relations 

SF = Fr{R,G)SR + Fg{R,G)5G, (27) 
d(^^R{R,G)SR + (G{R,G)SG. (28) 

In de Sitter space we have additionally 

G = AH^R, (29) 

which implies that 5G — AHq6R. Using the last relation 
together with (P7|) and (1^51) (and so limiting ourselves to 
de Sitter backgrounds) condition (|26|) becomes 

{/br + ^HlfRG + ^hUgb. + IGH^faa) SR = 0. (30) 

If f{R, G) is an analytic function we have /^g — fcR, 
and requiring that the above equation is valid for any 
variation SR (see the discussion on f{R) in Section llll Ap 
we arrive at 



The above equation is a second order partial differential 
equation with constant coefficients, for the class of func- 
tions / = f{R, G) that give a vanishing anisotropic stress 
in de Sitter space. Its general solution is 

f{R,G) = hm+Rf2m, (32) 

with n = R — G / {AHq), and /i, /2 arbitrary but analytic 
functions of f2. Models with vanishing anisotropic stress 
in de Sitter space need to be of this form locally near the 
de Sitter point. 

We specify the function /(-R, G) in the action, which 
is agnostic of quantities like Hq. For this reason it is 
preferable to consider a more general class of models with 

(33) 

with M a parameter with mass dimensions, so that 
jj(cff) _^ g corresponds to the special case of a model 
with a de Sitter expansion rate of i?o = M/2. As we 
will also discuss later on, the mass parameter M controls 
which of the two contributions in f{R, G) dominates. 

Assuming that the de Sitter point exists and is stable, 
we see that it is in principle possible to find a non-trivial 
class of /(-R, G) models that give exactly zero anisotropic 
stress in de Sitter space at all times, by selecting a model 
in the class ([5^ . However, as we will see by studying 
the stability of de Sitter space below, the case M — > 2i7o 
corresponds to a singularity for the actual model, and 
therefore the model cannot be viable. Furthermore, we 
will see that the anisotropic stress cannot become arbi- 
trarily small, since this will cause unstable behavior for 
the curvature perturbations. 

IV. ANISOTROPIC STRESS AND STABILITY 
FOR A DE SITTER BACKGROUND 

There are different stability criteria that a gravita- 
tional theory aiming to describe the late time acceler- 
ation should satisfy, each leading to a different condition 
for the form of the function f{R, G). At the background 
level, a viable model should give rise to sufficiently long 
radiation and matter e ras, as well as a transition to a 
stable de Sitter era |36l - [38| . Furthermore, avoidance of 
singularities and of rapid collapse of perturbations (posi- 
tivity of the sound speed) as well as agreement with local 
gravity constraints should be ensured [39l - l4l| . Of great 
importance is also the absence of ghost like degrees of 
freedom [42|-|4^. For the class of /(i?, G) models in a de 
Sitter background the latter requirement translates into 
/fl(i?,G) >0. 

Modified gravity models of the type f{R) ov R + /(G) 
suffer from a curvature singularity at very early times of 
the cosmological evolution p5l - |49| '*. The latter singular 



Irr + ^HlfRG + IQH^GG = 0. 



(31) 



^ Other types of singularities in tlie context of modified gravity 
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behavior can lead to oscillations of the scalar degree of 
freedom with infinite amplitude and frequency. As ex- 
plained in Ref. [i^, the singularity lies at a finite field 
value and energy level and therefore is easily accessible. 
We will see in the following that this singularity is a fea- 
ture of f{R, G) models as well. 

In this paper we are interested in the classical stability, 
and particularly its connection to the effective anisotropic 
stress. As we will show and discuss below, the attempt 
of turning off or making sufficiently small the effective 
anisotropic stress for a de Sitter background leads to se- 
rious stability problems that question the actual viability 
of models with vanishing n^'^*^^ . 



A. Existence of a de Sitter point 

Since we will specifically study the behaviour near the 
de Sitter point, it is necessary that this solution exists 
for the models of interest. De Sitter space is a vacuum, 
maximally symmetric space described by the conditions 

H = Ho = constant >0, R = G = F = i = 0. (34) 

Furthermore, in maximally symmetric spaces any curva- 
ture invariant can be expressed as a function of the Ricci 
scalar, and particularly for the Gauss-Bonnet term we 
get 



G 
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(35) 



We can derive the condition for the existence of the de 
Sitter point by taking the trace of the equations of motion 
([3]) and using relations (l34l) . ([SS]). to arrive at 



(36) 



F{R)R + 2G{R)^{R) - 2f{R) = 0, 



where everything is assumed to be expressed in terms of 
the Ricci scalar and evaluated on de Sitter space. The 
cases ^ = and F = 1, give the relevant conditions 
for f{R) and R + f{G) gravity respectively. Solving the 
algebraic equation given above, we get the de Sitter point 
solution, which in general is not unique. Minkowski space 
corresponds to the special case of Rq = Hq = 0- 

For the models of the type ([5^ , we find with the help 
of Eq. (pS)) that the de Sitter point is given by solutions 
of the equation 



fliu)+uf2{u)=0. 



(37) 



and R = 2u. The next step in our analysis will be the 
study of the stability of de Sitter space at both homoge- 
neous and inhomogeneous level. 



have been earlier observed in Ref. [50[ . For a discussion on 
various type of singularities in the same context see for example 
Ref. [sill and references therein. 



B. Homogeneous perturbations 

Now we turn to study the stability of the de Sitter so- 
lution, first with respect to homogeneous (background) 
perturbations. As we will see, there is a strong link be- 
tween effective anisotropic stress and stability in modified 
gravity models. 

Let us consider the t — t component of the Friedman 
equation ^ and perturb it linearly around the de Sitter 
solution H ^ Ho 



H{t) = Ho + SH{t). 



(38) 



Under perturbation 
reads as 



the perturbed function f{R, G) 



f = fo + Fo5R + ^oSG, 



(39) 



and similar expressions hold for the other quantities of 
interest. The explicit formulas and calculations for any 
space can be found in Appendix [K\ 

Now, evaluating relations (jA13p - (jA15|) and using con- 
ditions p4p . we can write the linearized perturbed mod- 
ified Friedman equation (jlip in the form 



CiSH + C2SH + C3SH = 0, 



(40) 



with the constants Ci, C2 and C3 defined in Appendix [Xl 
There is no constant term since we know that de Sitter, 
SH = 0, is a solution. This equation then admits an 
exponential solution of the form 



6H = A. 



with 



F 



and 



Fr + %Hl (Fg + 2HUg) 



(41) 



(42) 



(43) 



where we dropped the subscript "0" from Fr e.t.c for 
simplicity. 

From solution (|^^ . we can read off the condition for de 
Sitter space stability with respect to homogeneous per- 
turbations: 



F 



3 [Fr + 4ff2 {2Fg + 4i/2^G) 



AHl > 0. 



(44) 



The latter condition ensures that the de Sitter point is an 
attractor for the particular /(-R, G) model under study, 
which is important for the viability of a cosmological 
model of gravity. The limit ^ — >■ in (j^J) gives the 
corresponding condition for f{R) gravity, that has been 
derived before in Ref. [53 , 



F 
3F^ 



AHl > 0, 



(45) 
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while when _F — > 1 we get a similar condition for the 
R + f{G) models also derived in Rcf. ^48.] 



1 



- AH?: > 0. 



(46) 



The stability condition (H^ is general, but now we can 
check what it tells us for the class of models that give 
a zero anisotropic stress, described by equation (j33p as 
M'^ — )■ 477q. We can see that in this case necessarily 
cj — > 0, and so the eigenvalues ([^^ tend to infinity^. 
In particular, when uj is exactly zero, which corresponds 
to the case of a vanishing anisotropic stress, it is not 
possible to reach the de Sitter state without triggering 
a singularity in the model: the quantity Ci = ISHuj 
in Eq. goes to zero as we approach de Sitter, 

together with C2 — >• (see appendix In general 

SH ^ [Cj,/ Ci)5H — 00 which requires 5H to diverge 
in order to satisfy the evolution equation, except possi- 
bly for a lower dimensional and thus infinitely fine-tuned 
set of trajectories in specific models. We will show and 
discuss this explicitly in section |VI] considering examples 
for particular /(i?, G) models. 

Additionally, if the effective anisotropic stress is not 
exactly zero, but very small, then the rapid background 
oscillations lead to potentially large time derivatives of 
5H, which makes the linear analysis unreliable, i.e the 
evolution becomes non linear. We give an explicit exam- 
ple in section IVI Al 

Similarly, in the f{R) and R + /(G) cases, where the 
zero anisotropic stress condition was that /rr^R) = 
and /gg(G) = respectively, conditions pS)) and dH]), 
give the obvious result that one gets infinities when trying 
to suppress the extra degree of freedom. The difference 
with the more general f{R, G) models is that the singu- 
larity appears for a finite value of the mass parameter 
M of the model, while in f{R) and R + f{G) the same 
happens for rather trivial cases. We conclude therefore 
that a /(i?, G) type model that has no anisotropic stress 
in a de Sitter background cannot dynamically reach this 
background solution. 



C. Inhomogeneous perturbations 

In this subsection we will study the behavior of in- 
homogeneous perturbations in de Sitter space and we 
will first show that the stability condition coincides with 
the stability condition derived in Section [IV Bl on homo- 
geneous perturbations. This was found to be true for 
f{R) models in Ref. 



52| and in general is not true for 



scalar-tensor models. We will then make the relation be- 
tween anisotropic stress and stability clear by studying 



^ Since the no-ghost condition requires that F > 0, the question 
whether the background solution moves towards or away from 
de Sitter depends on whether uj — > 0"*" or a; — > 0^ . 



the evolution of the perturbations. The full set of per- 
turbation equations together with some useful relations 
can be found in the Appendix |B] 

We follow [23] and choose the gauge invariant expres- 
sion 



1 



^^ — [5F + ml6(\ 



(47) 



for the gravitational potential, as it reduces to (j) in the 
Newtonian gauge and remains well-defined for a de Sit- 
ter background. For that background, we find that the 
potential is just given by 



RR 



^Hq/rg 



IGH^fGc) SR 



2F 



(48) 



where we used the fact that Jrg — fcR and that in de 
Sitter space we have 5G = 4H^6R. From condition (|30|) 
we see that in de Sitter space and for models that have 
no anisotropic stress, 3> is necessarily zero. However, let 
us assume that we are not exactly in this limit. Then by 
substituting the expression of SR in terms of the gauge 
invariant $, relation (jBll|) . we arrive at the evolution 
equation. 



with 



'off 



n2 



F 
3w 



'off 



Am 



* = 0, 



(49) 



(50) 



with u! defined in (|43p . and a(t) cx exp(_ffof). m^g is the 
effective mass of the Klein-Gordon type equation for the 
scalar perturbation in de Sitter space, and has a purely 
geometrical origin. Equation p9)) reduces to that of f{R) 
and R + f{G) for the limits of ^ ^ and — >■ 1 respec- 
tively. 

As fc — > 0, the requirement for superhorizon stability 
dictates that the effective mass is positive. 



'eff 



> 



(51) 



which leads to the same stability condition as derived 
before with the homogeneous analysis, equation (|44l) . 
Therefore, the two stability criteria, with respect to 
homogeneous and inhomogeneous perturbations respec- 
tively, lead to the same conditions, as it is the case for 
f{R) gravity as well f52j . 

Turning back to the effective anisotropic stress, we can 
see that considering again the class of models found in 



and requiring M — >■ 2Ho (11 



(off) 



denominator of ([50)) go to zero so that 



0), will make the 



lim 

M^2Ha 



w^O \ 3aj 



= lim Am = ± 



ICX), 



(52) 



depending on the sign of w as it approaches zero. In 
the case of positive infinity the stability condition is not 
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violated, while the minus infinity will obviously violate 
the stability condition, as it would make the square of 
the effective mass negative (tachyonic). 

The effective mass going to infinity means that the 
scalar degree of freedom becomes frozen and so it is ef- 
fectively suppressed. This is also the case in the spe- 
cial cases of f{R) and R + f{G) gravity, as can be seen 
by inspection of equations (HSI) and for fuji — 
and foG respectively. However, here the singularity 
appears in a non trivial way, i.e for a critical value of 
the mass parameter M where the two different contribu- 
tions, i.e the R- and the G- contribution in (pS)) balance 
each other. By consequence, in f{R, G) type models, 
the anisotropic stress is related to the extra scalar de- 
gree of freedom of the theory. As the same happens in 
scalar-tensor models (e.g. equation (43) of [1^), and also 
in DGP where the absence of anisotropic stress requires 
the crossover scale to diverge, rc — > oo, which effectively 
restores GR, we conjecture that this is a quite general 
feature of modified gravity models. In addition, in the 
/(i?, G) case, turning the anisotropic stress off (or trying 
to make it sufficiently small) has a direct impact on the 
stability and time evolution of the model. 

To see what happens when the mass diverges, it is pos- 
sible to study the solution of the evolution equation (^5)) 
using a WKB approximation. We discuss the procedure 
in more detail in appendix [Cl where we show that the so- 
lution in this regime, and for a sufficiently large effective 
mass TOeff , is approximately given by 



$(t) « ^ C± exp 



-HqI ± irricst 



(53) 



with C± constants and Hq > 0. From the above solution 
it can be seen that the frequency of the oscillations is 
proportional to TOeff- Suppressing the anisotropic stress 
leads to a very large effective mass and thus to a very 
rapid oscillation of $. Although we have shown this here 
only for the de Sitter limit, we expect that the result 
is more general, and similar oscillations have been seen 
for example in Ref. (3^ during matter domination for 
numerically reconstructed f{R) models which mimic GR 
at early times. 

From relation ((50|) it can be seen that a sufficiently 
large effective mass corresponds to a sufficiently small 
anisotropic stress. However, the curvature perturbation 
6R (or SG) has an amplitude that is oc m^g. 



SR{t)=6{mls + 4H^)m 



(54) 



and as m^g 3> 1 one can get very large curvature pertur- 
bations. The latter behavior, occurring while we try to 
suppress the effective anisotropic stress (unless the ini- 
tial conditions are tuned appropriately), is very similar 
to the one caused by the singularity found in Starobin- 
sky's "disappearing cosmological constant" model, Refs 
|45| and [43. In that case, the singularity appeared in 
the high curvature limit of the particular model, while 
in our case it appears in the model space of different 



f{R, G) models respectively. The latter oscillatory be- 
havior endangers the stability of the actual model as has 
been pointed out in Refs [45] and [l^l , and for an explicit 
discussion on the subject the reader is referred to Refs. 

SSI- 

Another interesting aspect of the models of the type 
/(ri) concerns the sound speed. The propagation speed 
in de Sitter space equals the speed of lig ht (c^ — 1). 
However, using the formula derived in [23| we find that 
the sound speed in a general background is given by 



8H 



4ff 2 _ A//2 



^ (55) 



where 7 = jjp- is a dimensionless parameter (constant 
in a de Sitter background) .7^1 implies that the Ricci 
scalar part of the f{R, G) contribution to the anisotropic 
stress dominates, while for 7 ^ 1 the Gauss-Bonnet part 
is larger. 7 = 1 corresponds to the case where the two 
contributions in f{R, G) models become equal and can- 
cel. 

We can calculate H from the equations of motion, and 
for this particular class of models we get 



H 



(56) 



which can then be substituted in (|55p . However, con- 
sidering an expansion characterized by an effective WcS, 
equation ([22)) . the sound speed takes the form 



3(1 + Wcff) 
1-7 



(57) 



Assuming a background with Wcs 7^ — 1 , we immediately 
see that as 7 — >■ 1, — 00. The sound speed becomes 
negative for 7 < 1 (Gauss-Bonnet part dominates) and 
positive for 7 > 1 (Ricci scalar part dominates) respec- 
tively. The value 7=1, which corresponds to the effec- 
tive anisotropic stress becoming zero is the critical value 
where the sound speed diverges and changes sign. In 
other words, if one wishes to enforce Cg < 1 then one 
has to ensure that the model lies sufficiently far from 
the reg ime where the two contributions balance (see e.g. 
[53II54I for a discussion on possible issues of superluminal 
propagation of perturbations). 



V. GENERAL AND MATTER-DOMINATED 
BACKGROUND 

In this section we extend the analysis to a general back- 
ground evolution, and then consider specifically the im- 
portant case of matter domination. In general we have to 
consider Eq. ((23)) . In this equation, SF and are func- 
tions of SR and SG through Eqs. ^ and These 
in turn can be expressed in terms of the metric pertur- 
bations, (j) a-iid -0, see e.g. [5^. In the small-scale limit. 
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k aH, we find that (/> = tp implies 

fRR + 16iH^ + H)iH^ + 2H)fGG + 4.{2H^ + 3H)fRG - 0. 

(58) 

In order to re-transform this condition into one involving 
only R and G, we can eliminate H and H with the help 
of equations © and pUj) , 



12 



6 



(59) 
(60) 



Using this prescription we find for the general no- 
anisotropic-stress condition 



= fRR + ^{ -9G + 2R{R- VW^6G) ) faa 



R - 2 - 6G Jrg 



(61) 



While it is difficult to find general solutions, we can 
instead study the case for a background evolving with a 
given Woff, as defined in (1^ . We notice that in this case 
equation ([58|) can be written as 



= fnR-2H\h + 9w,s)fRG 

+ 8^4(2 + 9Weff(l+ii'cff))/, 



GG- 



(62) 



For Weff = — 1 (de Sitter expansion) we recover Eq. (1311) . 
while for WcS = (matter dominated expansion) we find 



Jrr - lOH^fRG + IQH^fGG = 0. 



(63) 



The Hubble parameter in the latter equation can be elim 
inated in favor of R and G using equations ([9|) and ^TU\ 
evaluated for a matter background, 



R = 3H^, G = -12H^, G = --R^. 



(64) 



We now try to construct an explicit example for a model 
that has no anisotropic stress during matter domination. 
For this purpose, we make an ansatz 



f{R,G) = R + G''(i(R) 



(65) 



Here we take /3 as an a-priori general function of R. In- 
serting this model into Eq. ()63p and using (|64|) we can 
re-express the condition in terms of R only. We find that 
P needs to satisfy the following differential equation: 



2n{n - 1)13 + bnR(3' 4- 2R^(3" = 



(66) 



This equation has clearly a power-law solution, j3{R) — 
cR"^, with 

mi,2 = i (2 - 5n ± + n{9n - 4)) , (67) 

and the general solution is of the form 

f{R, G)^R + ciG"i?™i + c2G"i?™^ (68) 



where = mi{n) is given by the equations for mi and 
1712 above. 

A successful model with zero anisotropic stress should 
at the same time satisfy the Friedmann equation as well. 
During matter domination we can write the latter as 







R^Irm 



'^G^/gg 



RGL 



RG 



RJr + Ig/g 
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If + lpoR. 



(69) 



Here we chose R and G so as to correspond to the partial 
derivatives, since the choice is not unique. The final term 
is due to pm{t) cx (x R. Inserting a model of the form 
(|68p but for a general exponent m, wc find the condition 

- 6to2 + m(7 - 6n) + [n - l)(12n - 1) = 0. (70) 

A model of this form that satisfies simultaneously ([S7)) 
and ([70)1 allows for a matter dominated evolution and 
contributes no anisotropic stress during that period. This 
is the case for 

n= — flliVUV TO= — f61±llV4T') (71) 
90 V y ' 180 V ^ ) \ ) 

where one needs to use either both positive or both neg- 
ative signs. An additional solution is given by m = and 
ri = 1, which is just GR. 

Therefore, there is at least one model in the context 
of /(i?, G) gravity that is able to give a zero effective 
anisotropic stress, in the subhorizon limit of a matter 
background. 

Let us now turn attention to homogeneous pertur- 
bations around the matter point, keeping the function 
f{R, G) in its general form for the start. In Appendix \K\ 
we calculate the evolution of homogeneous perturbations 
for a general expansion a{t) (x t^. For the matter case 
we get for p = 2/3 



with the effective mass defined as 

F_ ^ F 

3w ~ 3[Fr + 4F2 (2Fg + AH^^g)] ' 



cff 



(72) 



(73) 



Equation ([72]) can be solved approximately at the 
WKB regime using an iterative approach [4^ [43 , 



6H = 5H, 



(osc) 



(5^(md)- 



(74) 



(5ff(osc) is the solution describing oscillations of the scalar 
degree of freedom, obtained setting Sp^ = 0. (S^f(ind) 
denotes the matter induced part, which is obtained by 
turning off all the derivatives on the l.h.s of equation 
([7^ . We assume that \5Hi^asc)\ ^ '^^^(ind): so that the 
deviations from GR are sufficiently small. 

Stability in this case requires, apart from the no-ghost 
condition F > 0, that the effective mass is positive. 



TO 



off 



> 0. 



(75) 
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Let us turn attention to the oscihatory part of the so- 
lution ((7^ . It can be obtained using the WKB approx- 
imation, by assuming the solution is a slowly varying 
quantity in time. 



SH, 



(osc) 



(76) 



with ^ <C 1. Plugging above ansatz into ((72)) . and after 
some algebra, we find that 



6H, 



(osc) 



± cxp 



2 w 



X exp 



±i / dt' 



(77) 



with A± constants. Using the fact that H = Hi^{t) ~ 
2 /{it), and performing the integration in the first expo- 
nential we arrive at. 



5H, 



(osc) 



A± 



E (a;t3)l/2 



exp 



±i / dt' 



nics 



(78) 



The second integration can be performed after choos- 
ing a particular model. From (f78|) one can see that the 
amplitude of the oscillating solution grows as one goes 
backwards in time, which is exactly the behavior pointed 
out for f{R) models in Refs. |45l-l47j. and was due to a 
curvature singularity as explained in Ref. [l^]. There- 
fore, f{R, G) models suffer from the same problem too. 

For the model giving a zero effective anisotropic stress 
during matter domination, see relations (|68|) and (|71[) . 
one can check that both ui and m^g become complex, 
which is unphysical. What is more, the latter fact ren- 
ders the analysis of homogeneous stability for this model 
impossible, or at least highly non trivial. One could pos- 
sible seek to find a model without this peculiarity, how- 
ever we shall leave this for future work. 

We notice that the models of type (l32t that have no 
anisotropic stress during a de Sitter phase with specific 
expansion rate Hq = M/2 will pass through lo = and 
thus rrics — ?> oo in any background if the expansion rate 
H{t) crosses this critical value M/2. 

A different, general way to decrease the anisotropic 
stress is to move close to GR by decreasing the deviations 
from the extra f{R, G) contributions, which effectively 
implies 



/rr, Irg, Igg ^ 1- 



(79) 



In this case, we also make w small while F — > 1. Again 
this will lead to rapid oscillations, and we suspect that 
this is the reason for those seen in [3J] . Once the genie of 
extra degrees of freedom is out of the bottle, it is difficult 
to push it back in without further complications. 



VI. TOY MODELS 

In this section we will study the de Sitter behavior for 
some characteristic cases of the class of models found in 



321) . For the sake of generality we will consider 



G 



" = ^ + ^M^' 



(80) 



with e — ±1. The particular class of models with a van- 
ishing of the anisotropic stress, found in (1321) . correspond 
to e -1 and M 2Hq. 

First note that, for the class of models (|32|) . it is possi- 
ble to parametrize both the de Sitter existence and stabil- 
ity conditions in terms of the parameter 7, which controls 
the different regimes of the model. For simplicity and il- 
lustration let us assume that /2 = 0. Then, the de Sitter 
condition (1551) becomes 



/ 7 + e 



V27 + e 
with fn = /o(f^o), and 



/oOo - /(J7o) = 0, 



(81) 



(82) 



Furthermore, for the de Sitter stability condition (1FT|) we 
get 



7 



> 2 , 

7 + \27-1-e 



7 



(83) 



We will assume that fio > 0, 7 > and real. The limits 
7 — ^ 00 and 7 — >■ correspond to the pure f{R) and /(G) 
regimes respectively. 

In principle, we will assume that through (|81|) we can 
express Jlp in terms of 7 and the other possible param- 
eters of the model as f7o = fio (7, q), and then use ((83)) 
to get a constraining condition. 



A. f(Q.) = Q. + Q.\q{Q./c^) 

Here, c is a positive constant of mass dimensions. This 
model is able to re-produce a late-time acceleration, since 
at late times ^ 1, and the logarithmic term will dom- 
inate. In four dimensions the linear term 17 is essentially 
equivalent to the Ricci scalar R since the Gauss-Bonnet 
term does not contribute to the equations of motion. The 
absence of a Minkowski solution makes this model rather 
unrealistic. 

A non trivial de Sitter solution can be found using (|36p 



(84) 



and the Hubble parameter is then trivially given by ([82 
The stability condition ([55]) yields 



27^-1 
27 -I- e 



> 0. 



(85) 



For both branches, e = ±1, de Sitter space is stable when 
7 > V2/2. 
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To illustrate the singularity when trying to reach de 
Sitter, we set 7=1 and for simplicity c ~ -\/6e so that 
the de Sitter solution is given by ilo = 1- Expanding the 
equation of motion in we find to first order, 



2(l + 2(^ij)2) SH + 0{{5H)' 



0. 



(86) 



Only in the second order term a contribution SH{SH)^ 
appears. We notice that it is not possible to solve the 
first term for real SH, so that necessarily SH ex. 1/6H 
will diverge when we try to dynamically reach de Sitter. 
The only exception is 5H = 0, i.e. the solution that is 
always de Sitter. 

If we are looking at a model that is close to the critical 
case but has not quite zero anisotropic stress, 7 = 1 + e, 
then the stability analysis indicates that model is stable 
as long as e > —0.29, and indeed numerically we observe 
rapid oscillations of SH around de Sitter for most cases. 
For very small e we can get an idea of the model behavior 
by first solving the equation of motion for 6H and then 
linearizing the full equation with respect to SH, under 
the assumption that we will look at the evolution close 
to de Sitter. We can then expand the resulting equation 
in e, which to lowest order in e and for c = \/6e results 



(87) 



SH w —SH {2{SHy - 1 



We see that the second time derivative of SH will become 
very large as e — ^ 0, which shows again the presence of the 
singularity for the critical case of zero anisotropic stress. 
But we also see that even close to de Sitter the evolution 
of SH is governed by a non-linear differential equation 
since the time derivatives of SH will in general be large 
for models with small anisotropic stress and so cannot be 
neglected. A detailed study of how small the anisotropic 
stress can be made is therefore not straightforward and 
left to future work. 



B. 



f{Q,) = n + cn" 



In the context of /(i?) gravity, models of this type 
were sug gested as an explanation for late time accelera- 
tion j56j, |57[ with n < 0, while models with n > can 
lead to acceleration at early times and explain inflation. 
Furthermore, it was found that de Sitter space is unstable 
unless cn < [s^. Here, we assume that cn > 0, other- 
wise the no-ghost condition F > could be violated. 

The de Sitter point equation (|5T|) gives two solutions, 
namely fio = (for n > 0) which corresponds to 
Minkowski spacetime and a non trivial de Sitter one, 



7 



c(7(n-2) + e(?i-l)) 



l/(n-l) 



In order for 57o to be real and positive one has to ensure 
that the quantity in the denominator in the latter rela- 
tion is positive. We shall also require that the Hubble 



parameter, as given implicitly in relation (j82p . will be 
real and positive too. 

For both branches e = ±1, de Sitter is always unstable 
when n < 0. For n > 0, it is always unstable if e = 1, 
but for e = — 1, n > 2, the stability condition (|83| gives 



(89) 



with both f2 and Hq being real and positive. 

To avoid a superluminal sound speed, the model should 
lie in the f{R) regime, characterized by 7 > 1, which is 
satisfied here, as the right branch of above inequality 
approaches the value 1+ as n — )• 00. Further, for n > 2, 
Minkowski space is always stable. 

To consider the equation of motion close to de Sitter, 
we set e = —1, 7 = 1 and choose c = — 6*^"'^""^ for which 
Ho = I. We also assume that n ^ 1. We again expand 
in SH. The lowest order equation becomes now 



(1 -I- 2n(SHf^ SH + {{SH)^) = 0. 



(90) 



Again the second derivative of SH appears only at or- 
der (SH)^. This time we can in principle make the 
first order term vanish for rt < 0, which would allow 
to cross SH = with a finite second derivative. How- 
ever, there are two problems; Firstly, we can only cross, 
not move into and stay on SH ~ 0, since locally we need 
SH ^ (t — to)/\/—2n to avoid triggering the instabil- 
ity, and secondly this requires an infinite amount of fine- 
tuning in the initial conditions: we need to reach de Sitter 
at exactly the right speed, else we are either repelled, or 
a catastrophe engulfs the universe. So in reality again it 
is impossible to reach de Sitter dynamically. 



C. 



/(n) = f7-hcoA((l + ^)-"-l) 



This is a straightforward generalization of Starobin- 
sky's disappearing cosmological constant model [i^. It 
was proposed in the context of f{R) gravity as a late 
time acceleration model, that has a vanishing cosmolog- 
ical constant in Minkowski spacetime. It is trivial to 
check that Minkowski, /(O) = 0, is indeed a solution, but 
unstable since fnn{0) < 0. 

The model is characterized by three parameters, cq. A, 
7 > 0. 

From the de Sitter point equation, one can find an 
expression for A as a function of 7, and xi = Qq/cq 



A = 



xiigo - + xf) 



2\n+l 



[xl{2ngo + 1) - {1 + xl)^+^ + ly 



(91) 



where 170 = (7 + £)/(27 + e)- Taking the limit 7 — >■ cx) 
in the above expression one recovers the one given in 
Starobinsky's paper Ref. [45| . 

Let us assume that cq is of the order of the de Sitter 
scale, rio/co = xi = 1. The de Sitter stability condition 
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then reads 



2(7 + e)n^ + (27 + l)n + (27 + e)(l - 2") < 0, (92) 



For n = 1 de Sitter is stable if 

27 



27- 



>0, 



(93) 



which is never satisfied for both branches e = ±1. How- 
ever, choosing xi — 1/2, n = 1, we find that de Sitter is 
stable for e = 1 and 7 > 1/2, as well as for e = —1 and 
7 > 0.68. 

Stability can be established for a wide range of the 
model parameters, but that would require a detailed ex- 
ploration of the parameter space of {co, A, n}, and we are 
not interested in this here. 

For the critical case e = —1, 7 = 1, choosing the A of 
([?!]) and in addition cq = 6Hq for simplicity, we find to 
first order in SH 



[h^ - 2n{5H)' 

(2" - mt, 



SH + O {{5H'f) = 0. (94) 



This equation is of the same kind as the one found for the 
previous toy model, and it leads to the same behaviour. 
The special case n = 1 leads to the equation {SH)'^6H — 
0, which prohibits any crossing of SH = as otherwise 
SH has to diverge. 



VII. CONCLUSIONS 

In this paper we study the anisotropic stress in f{R, G) 
type modified gravity models. We investigated the pos- 
sibility of finding models that are able to mimic GR at 
least in the sense that they do not create an additional, 
effective contribution to the anisotropic stress, i.e 4> = 4' 
in the Newtonian gauge. For the needs of our analysis, 
we also derived the necessary background stability con- 
ditions. We started by considering the case of a de Sit- 
ter background, since this allowed us to find the general 
class of models with vanishing anisotropic stress. The de 
Sitter case is in addition interesting as current observa- 
tions indicate that the Universe is approaching this state. 
We further considered the general case in the small-scale 
limit, and in more detail the case of a matter dominated 
expansion. 

We find that for de Sitter expansion, the anisotropic 
stress is inextricably linked to the presence of an extra 
scalar degree of freedom. Suppressing the effective, geo- 
metric anisotropic stress is equivalent to suppressing the 
extra degree of freedom, which either requires the model 
to revert back to GR or else leads to an instability in 
the background evolution. In addition, it leads to prob- 
lematic effects like rapid oscillations of the gravitational 
potential and the curvature perturbation (with possible 



runaway production of scalar particles) and superlumi- 
nal propagation of the perturbations. The same prob- 
lems appear when one tries to generally decrease the ex- 
tra degrees of freedom through a model reconstruction, 
in order to obtain an evolution similar to GR. We think 
that this has been observed for numerically reconstructed 
f{R) models in a matter dominated background [s^, in- 
dicating that it is more general and not restricted to de 
Sitter. 

Furthermore, our stability analysis reveals that the 
curvature singularity present in f{R) models [45 4^ ap- 
pears in the more general f{R, G) case as well. What is 
more, its unwanted effect on the behavior of curvature 
perturbation is amplified for all models that try to sup- 
press the anisotropic stress by decreasing /flij, Jrg and 
fcG- In these cases we find rapid curvature oscillations 
with arbitrarily high amplitude as — t/i — > 0. 

In the case of a pure matter dominated background, we 
were able to construct an explicit model that gives a zero 
effective anisotropic stress in the subhorizon limit. At 
late times, when the gravity modifications are expected 
to appear and the evolution ceases to be matter domi- 
nated, this model will no longer give (j) — Tp. This could 
possibly be avoided by constructing such models for a 
whole expansion history including late-time accelerated 
expansion. However, such a procedure would necessarily 
involve significant fine-tuning as changes in the expan- 
sion rate would have to coincide with changes in the be- 
havior of the function f{R,G), which would in general 
depend sensitively on initial conditions. This appears to 
be rather difficult to construct. In addition, as discussed 
above, such a model would not be able to reach the de 
Sitter state without encountering a singularity. 

While the link between effective anisotropic stress and 
the scalar degree of freedom of the theory was studied 
here in the context of f{R, G) models, it is also present in 
scalar-tensor and DGP models: If a scalar-tensor model 
is coupled to the Ricci scalar in the action through F{(p)R 
then the anisotropic stress is proportional to {F' / F)6ip 
and the analogy to the f{R) case is obvious. In DGP, the 
effective anisotropic stress vanishes for Vc oc M^/M^ — >■ 
00 where M4 and Ms are the four- and five-dimensional 
Planck scales [H, [s^]- In this limit, the 5-dimensional 
part of the action is suppressed and only the usual 4D 
Einstein-Hilbert action remains. 

We conjecture that suppressing the effective 
anisotropic stress in modified gravity models is dif- 
ficult, if not impossible, to achieve in a realistic scenario. 
In models with a single extra degree of freedom that we 
looked at {f{R), f{G), scalar-tensor models and DGP) 
it is not possible at all to have no effective anisotropic 
stress except in the GR limit. In more complicated cases 
like f{R, G) it is possible to cancel the contributions to 
the effective anisotropic stress coming from several extra 
degrees of freedom, but this appears to be fine tuned and 
the resulting models tend to develop fatal singularities. 
This reinforces the role of the anisotropic stress as 
a key observable for current and future dark energy 
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surveys. While the observation of a strong anisotropic 
stress would point towards a modification of GR, the 
absence of anisotropic stress would present a significant 
challenge for modified gravity models and would require 
strong fine-tuning, which in turn favors scenarios where 
the dark energy is a cosmological constant or an extra 
minimally-coupled field with negative pressure. 
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Using relations given above, and after some algebra, 
the modified Friedman equation becomes 



with 




-Ci^G)SG+C^j^^SR+ 




Spi = 0, 


) 


= 6HF + 


3F + 36H^t 




(A7) 


C{K) 


= SH^Fr 


+ 3HFr + 12H^R 


(A8) 


C{G) 


= 3H^Fg 


+ 3HFg + 12H^G 


- \yiG). 


(A9) 




= 3HFr H 


h 12H''£,R, 




(AlO) 




= 3HFg - 


h 12H^^G- 




(All) 



Substituting for the perturbations of 5R, 5G and their 
derivatives we arrive at 



Appendix A: Homogeneous perturbations of f{R, G) 

In this section we will present the stability analysis of 
any fixed point of the the f{R,G) Friedmann equation, 
using homogeneous perturbations around the relevant so- 
lution. Our starting point is the t—t equation ([TT|). which 
for convenience we reproduce it here again. 



3H^F 



3HF 



1, 



(Al) 



li H = H{t) is a solution of above equation then per- 
turbing around it as H{t) — >• H{t) + SH{t), and keeping 
up to first order terms we get for the curvature scalars 
and their first time derivatives respectively 



R^ R + 6 [aH5H + 5H^ , 
G ^ G + 24 \2{2H^ + HH)5H + H'^SH 



(A2) 
(A3) 



The next step is to perturb the modified Friedman equa- 
tion pT|) . Particularly, the scalar potential becomes 

V + RSF + GS^ 

= V + {RFr + G^R,)SR + {RFg + G^g)SG 

= V + V^R,)SR + V^G)SG, (A4) 

and after evaluating the scalar field perturbations, it 
takes the form 

+ 6 (V(fl,) + 4i7V(G)) 6H + 2A{HV(R) 
+ 2{2H^ + HH)V(G))5H, (A5) 

where subscripts in brackets simply denote indices, while 
those outside brackets denote derivative with respect to 
the corresponding variable. 



Gi6H + C2SH + C3SH - Sp, = 0, 



with 

Gi = 6 
G2 = 6 

G3 = 6 



G, 



(R) 



(A12) 



(A13) 



C(R) + '^C\g)H^ + ^C^j^^H + 16G(^^(ij3 + HH) 

(A14)" 



1 ■ 



HF+-F + QH'i + AC(R,)H + 4G(^)i/ 
SC(G){'^H^ + HH) + 8C^^^{6H^H + + HH 



(A15) 



Defining uj = Fr + 8H'^{Fg + 2H'^^g), the generalisa- 
tion of Eq. (|43l) for arbitrary H, we find that always 



Gi = 18Huj. 



(A16) 



For a polynomial background expansion, described by 
a{t) (X tP, the other coefficients become 

= — [pio + 8H\1 + 3p){^R + 2H^^g) + (1 + '■^p)HFr] 

(A17) 



G3 = 3^^F + 2HF +12H 
4 



-H\iR + AHHg) - ^Hu 
p J J 



(A18) 



For a de Sitter expansion, a{t) oc exp[HQt], and H 
Hq = const., we get 



G2 — SHqCi 
F 



3w 



- 4i?o' Gi 



(A19) 
(A20) 
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Appendix B: Inhomogeneous perturbations and de 
Sitter stability 

The general metric elenient for scalar perturbations 
around a flat FRW background reads 

ds^ = - (1 + 2a)dt^ - 2a{t)d,P dt dx' 

+ a{tf {5,j - 2(t)5,j + 2d,dj-i) dx'dx^ . (Bl) 

The general form of scalar perturbation equations 
around FRW for f{R, G) models can be found in Ref. 
[27| . Here, we shall present the full set of equations for 
the case of de Sitter space only. 

Before we proceed, let us define the gauge invariant 
variable $ as 



5F + m^5^ 

2F ' 



(B2) 



with H = Hq as well as the rest of the background quan- 
tities evaluated on the de Sitter point. The perturbation 
equations then read as 



-3iJ(/) = 3i/$+(^-3ij2)$, (B3) 



Now, using equation (IB4p in (jB3[) we get 

$ = + i/x, (B9) 
which can be inserted into (jBSP to give 



5- = 0. 



(BIO) 



Using equations (jB4[) . (|B5I) as well as (jB9|) we can re- 
express the curvature perturbation in terms of the gauge 
invariant potential $ 



6R=-Q 



$ + 3i?$ 



(Bll) 



Appendix C: Sub-horizon solution for $ in the WKB 
approximation 

Considering the evolution equation (|49|) in de Sitter 
space for the gauge invariant potential $, we assume a 
solution of the form 



(CI) 



X + iJx + - V = 2$, 



6R = -2 



12i7> 3<?!) -f l2H(j) + 3i/V 
^ (X + 2Hx + 20 - V) 



5G = -8 



12H4a - ZH^ij) + 3ij3d - 12ij3<; 
+ {2Hx + X - 2</> - a) 



(B4) 



(B5) 



(B6) 



(B7) 



Equations (|B3|) . ((B4)) and ((B5|) correspond to the 00, the 
Oi and the 7^ j) components respectively. Particu- 
larly, equation (jB5j) is the anisotropy equation, and the 
choice of variable $ is now evident: it is the r.h.s of the 
latter equation, describing the effective anisotropic stress 
in de Sitter space, $ — 11*^'^^^ , and therefore is gauge in- 
variant. 

In order to re-express above equations in terms of 
gauge invariant variables only, we need a second gauge 
invariant variable apart from $. Following [27j we define 



with C a constant, and 9{t) <Si 1. Then, we can calculate 
that 



$(t) « ^C±exp i / dt'9^{t') 
= ^C±exp|y dt' 



(C2) 



with 



3Hn 



A^^, B{t)^-^A (fc2e-2«o* + ,^2^) _ ml (C3) 



From solution (|C2I) we can calculate the limit when 
m^g ^ 1, which is the case when 11^''^' — > 0. In this case 
we have, 

B(t)«TOeff, B(i)»0, (C4) 
and the solution is approximately given by 



^{t) « 5] C± exp 



(C5) 



* = $ + - i/x- 



(B8) 
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